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PACS 67 . 85 . -d - Ultracold gases, trapped gases 

PACS 03.65.Vf - Phases: geometric; dynamic or topological 

PACS 73 . 43 . -f - Quantum Hall effects 

Abstract -We describe a simple scheme by which "optical flux lattices" can be implemented in 
ultracold atomic gases using two-photon dressed states. This scheme can be applied, for example, 
to the ground state hyperfine levels of commonly used atomic species. The resulting flux lattices 
simulate a magnetic field with high mean flux density, and have low energy bands analogous to 
the lowest Landau level. We show that in practical cases the atomic motion significantly deviates 
from the adiabatic following of one dressed state, and that this can lead to significant interactions 
even for fermions occupying a single band. Our scheme allows experiments on cold atomic gases 
to explore strong correlation phenomena related to the fractional quantum Hall effect, both for 
fermions and bosons. 
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There has long been an endeavour in the field of ultra- 
cold atomic gases to achieve an experimental regime in 
which strong correlation phenomena associated with the 
Pactional quantum Hall (FQH) effect can appear [T}{3]. 
Central to this goal is the formation of Landau levels: 
the highly degenerate energy levels of a two dimensional 
charged particle in a uniform magnetic field. Owing to 
this degeneracy, particles occupying a single Landau level 
are susceptible to interparticle interactions. These can 
give rise to strongly correlated FQH states when the par- 
ticle density is comparable to the density of magnetic flux 
quanta n^. Existing experimental studies of FQH states 
in semiconductors involve fermions (electrons). However 
it is expected that related strong-correlation phenomena 
can arise also for bosonic species of cold atoms [3]. 

The formation of Landau levels for neutral atoms re- 
quires the creation of an effective magnetic field. To date, 
this has been achieved in experiments on ultracold atoms 
cither by using rotation [HGI|, or by using optical dress- 
ing to generate effective gauge fields [3]. However, in both 
cases the magnetic flux densities achieved so far are 
too low to bring a large atom cloud into the fractional 
quantum Hall regime. Theoretical proposals for methods 
to generate effective magnetic fields for atoms on deep 
optical lattices hold promise to achieve high flux densi- 
ties, on the order of the inverse optical wavelength squared 
~ 1/A 2 (for reviews of these proposals see, e.g. [5JE]). 
At these flux densities, interaction energy scales would be 



sufficiently large to allow experimental studies even in the 
FQH regime. 

A recent theoretical proposal has shown that effective 
magnetic fields with high flux density ~ 1/A 2 can also 
be generated for optically dressed states without the use of 
deep optical lattices. In this so-called "optical flux lattice" 
scheme [7] a periodic effective magnetic field with non-zero 
average is generated by arranging that the spatial depen- 
dence of the Rabi coupling and a state-dependent poten- 
tial are in register with each other. The implementation 
proposed in Ref. [7J involves a single-photon coupling, suit- 
able for Ytterbium or alkaline earth atoms, and requires 
locking the positions of standing waves of a laser at the 
coupling frequency with those of a laser at an "anti-magic" 
frequency. 

In this paper, we show how optical flux lattices can be 
generated by using two-photon dressed states. Our scheme 
also requires two laser frequencies, but it is more robust 
than that of [7J: the Rabi coupling and state-dependent 
potential are automatically in register. Importantly, our 
scheme can be applied to commonly used atomic species, 
including alkali atoms in hyperfine levels of any angular 
momentum F . We describe the bandstructurcs for rep- 
resentative cases F = 1/2 and F = 1. We show that 
there appear low energy bands that are analogous to the 
lowest Landau level: of narrow width in energy and with 
nonzero Chern number Owing to the narrow energy 
width, particles occupying these bands are susceptible to 
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interactions and to the formation of strongly correlated 
FQH states. We show that, even for fermions interacting 
with contact interactions, there remain significant inter- 
particle interactions within this low energy band. Thus, 
our scheme will allow experiments on cold atomic gases to 
explore strong correlation phenomena related to the frac- 
tional quantum Hall effect for both fermions and bosons. 

In the first part of this paper we consider an atomic 
species with a ground level g of angular momentum J g = 
1/2. Examples of atoms in this category that have already 
been laser-cooled are m Yb or 199 Hg (level 6%) HQS]. 
The atoms are irradiated by laser waves of frequency ujl 
that connect g to an excited state e also with angular mo- 
mentum J e = 1/2. For Ytterbium and Mercury atoms, 
we can choose e to be the first excited level 6 3 Po entering 
in the so-called 'optical clock' transition. The very long 
lifetime of e (~ 10 s for Yb [TT] and ~ 1 s for Hg [TJ]) guar- 
antees that heating due to random spontaneous emissions 
of photons is negligible on the time scale of an experi- 
ment. Another possible choice could be 6 Li atoms, but we 
estimated in this photon scattering rate that is too 

large to maintain the gas at the required low temperature. 

We assume that the atomic motion is restricted to the 
xy plane and described by the Hamiltonian 
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2M 



1 + V(r) 



(1) 



where M is the atomic mass and p its momentum. The 
matrix V acts in the Hilbert space describing the internal 
atomic dynamics. For an off-resonant excitation, we can 
assume that the population of e is negligible at all times, 
so that V is a 2 x 2 matrix acting the g± manifold [13 . Its 
coefficients depend on the local value of the laser electric 
field, which we characterize by the Rabi frequencies n m , 
m = 0,±1, where mh is the angular momentum along z 
gained by the atom when it absorbs a photon. 

In order to increase our control on the spatial variations 
of V, we suppose that a magnetic field parallel to the z-axis 
lifts the degeneracy between the states g±. The resulting 
splitting S is supposed to be much larger than the K m 's. 
Hence for a monochromatic laser excitation at frequency 

u>l, the off-diagonal matrix elements V-\ and V-+ are 

negligible compared to the diagonal ones. However we 
also assume that another laser field at frequency uj, + 5, 
propagating along the z axis with er_ polarization (i.e. 
m = — 1 with the notation above), is shone on the atoms. 
The association of this field with the 7r component (m = 0) 
of the light at wl provides the desired resonant Raman 
coupling between | g±) (figure [T^,). Using standard angular 
momentum algebra we find in the |<7-)} basis: 



V 



nK tot 

3A 



1 



h 

3A 



EkZ 



Ek 



«+ 



(2) 



Here n 2 ot = J2 m l K m| 2 , A = uj l - u> A , where u> A is 
the atomic resonance frequency, and we assume |A| 3> 
\5\, \K m \- The quantity E characterizes the field of the ad- 
ditional laser at ujl + 5. This beam is assumed to be a 




Fig. 1: (Colour on line) (a) A ground level with angular mo- 
mentum J g = 1/2 is coupled to an excited level also with an- 
gular momentum J e = 1/2 by laser beams at frequency ujl and 
ujl + 8. The Zeeman splitting between the two ground states 
g± is 5. (b) Three linearly polarized beams at frequency ujl 
with equal intensity and with wave vectors at an angle of 2n/3 
propagate in the xy plane. The fourth, circularly polarized 
beam at frequency ujl + 5 propagates along the z axis. 



plane wave propagating along z, so that E is a uniform, 
adjustable coupling. The a.c. Starkshift due to this addi- 
tional laser is incorporated in the definition of 5. 

We consider the laser configuration represented in 
Fig.fljb). The laser field at frequency wi is formed by 
the superposition of three plane travelling waves of equal 
intensity with wavevectors ki in the xy plane. We fo- 
cus on a situation of triangular symmetry, in which the 
three beams make an angle of 27r/3 with each other, 
fci = -k/2 (V3, 1,0), fc 2 = k/2 (\/3,-l,0) and k 3 = 
k (0, 1,0). Each beam is linearly polarized at an angle 9 
to the z-axis, which leads to 
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cos z + sin 8 (z x ki) 



(3) 



where k is the Rabi frequency of a single beam. In the 
following we denote V = Hk 2 /(3A) the energy associated 
with the atom-light interaction and e = E/k the relative 
amplitude of the ujl + S field with respect to the ojl field. 
The recoil energy Eji = h 2 k 2 /2M sets the characteristic 
energy scale of the problem. 

The coupling V is written in Eq. @ as the sum of 
the scalar part hK 2 ot /('SA) 1 and a zero-trace component 
that can be cast in the form W = a ■ B/2, where the 
&i are the Pauli matrices (i = x,y,z). For J? / and 
sin 2(9 7^ 0, the coupling B is everywhere non-zero. Sup- 
pose that the atom is prepared in the local eigenstate 
|x(^)) of W, with a maximal angular momentum pro- 
jection along n = —B/\B\. Suppose also that it moves 
sufficiently slowly to follow adiabatically this eigenstate, 
which is valid when V 3> _Er. This leads to the Berry's- 
phase-related gauge potential ih(x\^x)> representing a 
non-zero effective magnetic flux density |14j . For most 
optical lattice configurations, the periodic variation of the 
atom-laser coupling leads to a zero net flux of the effective 
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field through a unit cell of the lattice. Indeed the peri- 
odicity of the eigenstate \x) entails that the line integral 
of (x|V%) on the edges of the unit cell vanishes. However 
this reasoning may fail if the vector potential has singular- 
ities inside the cell: this is the so-called flux lattice scheme 

The singularities of (x\^X) occur at places where the 
g+-g~ coupling [oc kq, see Eq. ([2])] vanishes. Here, the 
zeroes of k$ are located on a honeycomb lattice, with a 
distance of 2A/(3v / 3) between nearest neighbours, where 
A = 27r/fc [TS]. The unit cell of this hexagonal lattice has 
an area A ce n = 2A 2 /(3v / 3), and contains two inequiva- 
lent sites that correspond to vortices and antivortices of 
(n Xl n y ), respectively. We will introduce later a gauge 
transform that allows one to use the same unit cell for the 
periodic coupling V. At the location of the (anti)vortices 
of (n x ,n y ), the modulus of B z cx |«+| 2 — |^-| 2 is maximal 
and its sign alternates, leading to n z = ±1. Because the 
sign of n z differs for vortices and antivortices, each of these 
regions wraps the Bloch sphere in the same sense provid- 
ing the desired rectification of the magnetic flux. The net 
result is that n wraps the Bloch sphere once within the 
unit cell, corresponding to = 1 flux quanta. This net 
flux Nff, = 1 is a robust feature of the optical coupling, 
not relying on fine-tuning of parameters or on the control 
of the relative phase of the laser beams. Reversal of the 
magnetic flux (N^ = —1) can be achieved by changing 
the sign of the detuning A, or replacing the cr_ laser at 
frequency u>l + 8 by a er + laser at frequency lul — 8. 

Our laser configuration is reminiscent of the scheme 
studied in |16j . However in that work only the three beams 
at frequency u>l and propagating in the xy plane were 
considered. The off-diagonal coupling in V is in this case 
V-\ cx k+Kq + KoK*_, and one can check the desired rec- 
tification of the magnetic flux at the zeroes of V-\ docs 

not occur. Therefore, although the Berry's-phase-related 
magnetic field found in |16j is locally non zero, the flux 
per unit cell vanishes. 

In the adiabatic limit, valid for V ^> En [5], the atoms 
are strongly confined to the sites where SK 2 ot /(3A) — | _B | / 2 
is minimum. Depending on the sign of A and the values 
of 9 and E, these sites coincide either with the vertices of 
the hexagons of the honeycomb lattice introduced above, 
or with the centre of these hexagons, forming thus a tri- 
angular lattice. Both of these tight-binding limits are sit- 
uations in which the number of flux quanta per hexagonal 
(rcsp. triangular) plaquctte is integer (rcsp. half-integer), 
i. e. time-reversal symmetry is preserved^. 

In particular in this adiabatic limit the various energy 
bands cannot have a non-zero Chern number [S]. However, 
an important aspect of the physics of optical flux lattices 
is that they do not rely on the adiabatic limit, and can 
be applied in regimes of intermediate coupling V ~ £r. 
Indeed the optical flux lattices described here, and those 



1 ln both cases, the unit cell of the lattice contains = 1 flux 
quantum. 



of [7], lead to bands with non-zero Chern number in the 
weak coupling limit V < E^. In this regime, a tight- 
binding model is inappropriate; departures from this (with 
next nearest neighbour tunnelling and beyond) allow the 
atoms to perform loops that enclose flux that breaks time- 
reversal symmetry. 

We have calculated the band structures for arbitrary 
V / Er by numerical diagonalization of the Hamiltonian 
|T|). To do so, it is helpful to expose the full transla- 
tional symmetry of the system. The matrix V is invariant 
under translations by the vectors A/\/3 (±1, y/3, 0), corre- 
sponding to a 'naive' unit cell of area 2A 2 /V3. One can 
reduce the size of the unit cell by taking advantage of a 
unitary transformation by U = exp(— ik^ ■ ra z /2). This 
state-dependent gauge transformation gives a transformed 
Hamiltonian 

H' = &HU = iP - dzhk " /2? + V', (4) 

2M ' v ' 

where V' = U'VU is of the same form as V, but with 
no replaced by k' = e~ 3 ' r «;o- The advantage is that 
the coupling V' involves only the momentum transfers 
K± t 2 = fci,2~ &3- Using the momentum transfers to define 
the reciprocal lattice vectors leads to the largest possible 
Brillouin zone, and smallest real space unit cell, which co- 
incides with the unit cell of the scalar potential described 
above, with area ^4 C cii = 2A 2 /(3-\/3). 

For the special case e = 9 = 0, the atoms only ex- 
perience the scalar potential oc |«o| 2 . For A > 0, this 
situation was investigated in detail in |15j . in connection 
with graphene physics. As explained above, |ko| 2 has min- 
ima at the sites of a honeycomb lattice. For V > -Er one 
therefore expects the atoms in both g± states to have the 
low-energy bands of the honeycomb lattice with two Dirac 
points in the Brillouin zone. One difference with respect to 
graphene is that, due to the state-dependent gauge trans- 
formation U, the momenta of the two spin states arc offset 
by ±k$/2. This has the effect that there are three special 
points in the Brillouin zone: a Dirac point for g + , a Dirac 
point for and a (coincident) Dirac point for both g±. 
The bands are shown in Fig. [5] with dotted lines for a cut 
through the Brillouin zone that passes through these three 
points. 

Perturbations to this decoupled limit cause the Dirac 
points to split, with consequences that depend on the 
symmetries that are broken and the relative sizes of 
the perturbations |T7]. For small non-zero 9, the terms 
|k + | 2 — |k-| 2 cx 9 2 in break the sublattice symmetry of 
the honeycomb lattice for both g± . Gaps then open in the 
spectrum at the locations of the Dirac points. For e = 
the lower two bands separate from the upper two bands in 
such a way that both pairs of bands are topologically triv- 
ial, that is each pair has net Chern number of zero. When 
both e and 9 are non-zero, the optical dressing leads to 
a net flux through the unit cell, indicating time-reversal 
symmetry breaking. Indeed, we find that for e, 9 ^ the 
bands can acquire non-zero Chern numbers. Specifically, 



p-3 



N.R. Cooper and J. Dalibard 




Fig. 2: (Colour online) Left-hand panel: Bandstructure for 
F = 1/2 for V = 1.8 E n along a path fc = k(K 2 - 2Ki) - fc 3 /2 
through the Brillouin zone. For e = 9 = (dotted blue line), 
the decoupled m — ±1/2 states each have two Dirac points. 
For weak coupling e = 6 = 0.1 (dashed green line) the Dirac 
points split in a manner that breaks time-reversal symmetry, 
giving the lower two bands a net Chern number of 1. For 
intermediate coupling e = 0.4, 9 = 0.3 (solid black line) the 
lowest energy band has Chern number 1. Right-hand panel: 
The density of states for V = 1.8 E R , e = 0.4 and 9 = 0.3. 



in the perturbative limit (e, 9 <C 1) the Dirac points split 
in such a way that the lower two bands have a net Chern 
number of 1 provided 9 2 /e is sufficiently small. When 
9 2 /e exceeds a certain value (~ 0.19 for V = 1.8 £/r) there 
is a transition to the topologically trivial case described 
above. Beyond the perturbative limit, as the couplings e 
and 9 increase, the splitting between the lower two bands 
increases and the lowest energy band evolves into a nar- 
row band with Chern number 1 . An example is shown by 
the solid lines and density of states in Fig.[2J for which the 
lowest band has a width AE ~ 0.1 The optical dress- 
ing ([2]) affords a great deal of freedom to tune parameters 
to reduce the ratio of the bandwidth of the lowest band 
to the gap to the next band. For example, for V = 2 E-&, 
9 = 7r/4, e = 1.3, the lowest band has a width of only 
AE ~ O.OIEr and is separated from the next band by 
about 0.4 E K (See Fig.©. 

It is important to emphasize that the formation of this 
narrow low energy band is not simply due to compression 
into a tight-binding bancj^. Rather it is closely related 
to the formation of Landau levels in a uniform magnetic 
field. A continuum Landau level is highly degenerate, with 
degeneracy equal to the number of flux quanta piercing the 
plane. Thus, for the flux density here, of one flux quantum 
per unit cell, a Landau level would have one state per 
unit cell: that is one band within the Brillouin zone. The 
lowest band of Figs. with its narrow width and Chern 
number of 1, is the optical flux lattice equivalent of the 



2 For a tight binding band in the limit of vanishing tunnel coupling 
when the Wannier states become exponentially localized, the Chern 
number of the band must be zero I18| , or a set of bands with net 
Chern number zero must become degenerate. 
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Fig. 3: Density of states for F = 1/2 for V = 2E K , 9 = tt/4, 
e = 1.3. The lowest band has Chern number 1, a width of about 
O.OIEr, and is well separated in energy from the next band. 
The density of states for the lowest band has been rescaled by 
1/10. 



lowest Landau level. 

The above scheme can be generalized to atoms of the 
alkali-metal family, whose ground state n Si / 2 is split in 
two hypcrfine levels I ± 1/2, where / is the nuclear spin. 
The laser excitation is tuned in this case around the res- 
onance lines D\ (coupling to nPx/2 with detuning Ai) 
and £>2 (coupling to nP 3 / 2 with detuning A 2 ). Let us 
focus here on the lowest hyperfine level F = I — 1/2. For 
the configuration of Fig.[T] the atom-laser coupling can be 
written 
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where A 1 



(5) 

F is the angular 



(l/3)Ar 1 + (2/3)A 2 -\ 
momentum operator in the ground state manifold in units 
of h, and 



B x + iB y — £ Ekq, 



B, 



M 2 



(6) 



with £ = (A^ 1 - Ai 1 )h/[3(F + 1)]. Under the unitary 
transformation U = cxp(— ik 3 ■ r F z ) the Hamiltonian 
takes a similar form to (U]), now with a z /2 replaced by 
F z , and again with a coupling V' in which kq is replaced 
by k' q = e~ lk3 ' T Ko giving the unit cell of the honeycomb 
lattice as before. Adiabatic motion of the atom still leads 
to a dressed state with angular momentum along the vec- 
tor n that wraps the Bloch sphere once in the unit cell. 
However, now the Berry phase acquired is larger by a fac- 
tor of 2F HH ■ Therefore, the unit cell contains N$ = 2F 
flux quanta. This increase of leads to an important 
new feature: a continuum Landau level now corresponds 
to N^, = 2F states per unit cell. Thus, the analogue of the 
lowest Landau level is a set of 2F low-energy bands with 
a net Chern number of 1. Spatial variations in the scalar 
potential and flux density will cause these bands to split 
and to acquire non-zero widths. 

We shall illustrate the physics for F > 1/2 by describ- 
ing the properties for bosonic atoms with F = 1. This 
is a very important case, as it applies for example to the 
hyperfine ground states of 23 Na, 39 K, and 87 Rb. In order 
to control the respective role of the various Zceman sub- 
levels, we also include in the Hamiltonian the additional 
Zeeman coupling to an external magnetic field 



(7) 
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Fig. 4: Density of states for F = 1. (a) V = 8£ R ,e = 0.8, 6> = 

0. 2, 71 = 0, 72 = —0.15 -Er. The two lowest bands have Chern 
numbers of —1 and 2. Taken together, these are analogous to 
the lowest Landau level (LLL). (b) V = 2 E R , e = 3, 9 = 0.5, 
71 = 10 £r, 72 = — 10 -Er. The system behaves as a two-level 
system formed from m — 0,-1. The lowest energy band is 
analogous to the LLL. 

where 71 and 72 correspond to the linear and quadratic 
Zeeman effects, respectively. For simplicity we assume 
that I Ai| <C IA2I and we take into account only the cou- 
pling between the ground state and the P1/2 excited state 
(Di line). As above, we characterize the atom- laser cou- 
pling by the angle 9 and the two dimcnsionlcss parameters 
e = E/k and V/E R , where V = ^k 2 /(3Ai). 

For a state-independent potential (e = 9 = 71 = 72 = 0) 
and positive detuning Ai, each of the three states \F = 

1, m) with m = 0,±1 experiences a honeycomb lattice 
potential and has two Dirac points. Now the offsets in 
wavevector are such that there are three special fc-points, 
at each of which two Dirac points coincide. Again the 
Dirac points are split by perturbations that break time- 
reversal symmetry (e, 9^0) and the lowest three bands 
can acquire non-zero Chern number. Differences from F = 
1/2 emerge in the regime of intermediate coupling. In 
particular we find a domain in parameter space where the 
two lowest energy bands are closely spaced, with narrow 
total width and a net Chern number of 1 [see Fig.HJa)]. 
This confirms the expectation mentioned above that the 
set of the 2F = 2 lower bands constitutes in the case the 
analogue to the lowest Landau level. One can also make 
use of the quadratic Zeeman splitting 72 =/= to arrange 
situations in which only two of the m-states contribute. 
This two-level system is then of the same form as F = 
1/2. In this case, parameters can be chosen to form a 
single narrow lowest energy band with Chern number 1, 
which by itself is analogous to the lowest Landau level. An 
example is presented in Fig.^b), showing similar features 
toF = l/2 (Fig.EJ. 

These narrow bands with non-zero Chern number are 
excellent systems with which to explore strong-correlation 
phenomena related to fractional quantum Hall physics [3J , 
both for fermions and bosons. Given the contact nature 
of interactions in ultracold atomic gases, one might think 
that the Pauli principle would prevent intcrparticlc inter- 



actions between fermions within a single band. Interac- 
tions would indeed vanish if the atoms were to move adia- 
batically, since then, at each point in space, only a single 
dressed state would be relevant. However, non-adiabatic 
effects lead to residual interactions between fermions. In 
a semi-classical picture, one can imagine two atoms that 
approach the same point in space from different directions. 
Since the motion is non-adiabatic, when they reach that 
point they can be in different internal states, so there is 
some probability that they coincide and interact. 

To understand quantitatively how non-adiabatic motion 
of the fermions leads to interparticle interactions within 
a single band, consider a system of spin- 1/2 fermions 
(m = ±1/2) with spin-independent contact potential 
V mi m> (r, r') = gw5 2 (r — r'). The mean interaction en- 
ergy for any Fock state is 

(Hi) = -92D /X/ [P-m,mRm',m' -\Pm,m'\ 2 ]d 2 r 

rn,m' 

= \ 9 2D J p 2 (r)[l-P 2 (r)]d 2 r (8) 

where p m ,m'(f) = ( c m( r ) c m'(r)) and (J m (r) creates a par- 
ticle at r in spin state m. In Eq. © we have introduced the 
local total density p and local polarization P (in units of 
h/2). In the adiabatic limit, the local polarization is every- 
where maximal, P = 1, so the interaction energy vanishes. 
However, for non-adiabatic motion the local polarization 
is reduced and the interaction energy is non-zero. 

We quantify the size of the interactions by considering a 
filled band. Then p m , m >(r) = £ fc [^m( r )Y ^m'( r ) where 
4>m(r) is the space- and spin-dependent wavefunction of 
the Bloch state k in this band. We define an interaction 
parameter I = (Hi)/Eq with E = (l/4)p2D / p 2 d 2 r the 
mean energy for P = 0. We have computed the total in- 
teraction energy for the lowest energy band in Fig. [2] for 
each of the three case^l of (e, 6). For e = 9 = the two 
spin states are decoupled, so the lowest band is unpolar- 
ized (P = 0) and 7=1. For weak coupling e = 9 = 0.1 
we find I — 0.95, showing a (small) suppression of interac- 
tions. With increasing coupling, the interactions continue 
to be suppressed. However, even for intermediate coupling 
(e = 0.4, 9 = 0.3), where the lowest band has become nar- 
row, the interaction parameter is I — 0.42, showing that 
the motion is still far from adiabatic. The ultra-narrow 
band shown in Fig. [3] has a small value for the interaction 
parameter I = 0.046. 

Since the low energy bands described above are anal- 
ogous to the lowest Landau level, leading candidates for 
strongly correlated phases include the Laughlin states, at 
mean particle densities of p = (1/3)N^/A ce ]i for fermions 
and p = (l/2)N < p/A ce n for bosons, with = 2F; or, if 
the Zeeman effect is used to isolate an effective two-level 
system, with = 1. Strongly correlated phases of this 

''Although in the first two cases there is an energy overlap with 
the second band, one can still compute the interaction energy for 
the Fock state in which the lowest energy band is filled. 
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kind are energetically competitive when the mean interac- 
tion energy is on the scale of the bandwidth AE ~ 0.1 Er. 
For bosons, taking £\ = g 2T) J p 2 d 2 r w (1/4)0 2 d-W|/A cc h 
as a typical scale for the interaction energy, we obtain us- 
ing the cell area given above E-J AE w (l/4)N^g, where 
g = Mgzo/fr 2 is the dimcnsionless interaction strength 
in 2D. With cold atomic gases confined in a strong one- 
dimensional optical lattice along the z direction, values of 
g can approach unity |19j . so that reaching the strongly 
correlated regime can be realistically envisaged. For spin 
1/2 fcrmions, the partial suppression of s-wave interac- 
tions discussed above somewhat complicates the realiza- 
tion of the correlated regime. However one can rely on 
Feshbach resonances to produce a strongly interacting 2D 
Fermi gas [50], and thus also in this case reach the point 
where (Hj) ~ AE. Note that in the case of 171 Yb or 
199 Hg, one has to turn to optically induced Feshbach reso- 
nances [2TJ[22] , because the spin of the ground state has a 
purely nuclear origin hence 'standard' magnetic Feshbach 
resonances do not occur. 

The appearance of fractional quantum Hall states is pre- 
dicted to have to various consequences for the ground state 
and excitation spectrum of the cold atomic g dis- 
cussed in the literature [THS]. Perhaps the most evident 
feature is a plateau (or plateaus) in in situ images of the 
density of a trapped gas, arising from the incomprcssibility 
of the FQH states. For the Laughlin states, these plateaus 
are expected at mean densities for which there is no com- 
peting Mott insulator state. 

In summary we have presented a simple scheme by 
which optical flux lattices can be implemented using two- 
photon dressed states. The method applies to com- 
monly studied atomic species, requires few lasers, and 
uses only well-established techniques. Our studies show 
that for practical parameters, the non-adiabatic charac- 
ter of atomic motion can lead to significant interactions 
even for fcrmions occupying a single band. This scheme 
should allow experiments on cold atomic gases to explore 
strong correlation phenomena in regimes of flux density 
previously unattainable. 
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